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Abstract 

The chromo-dielectric soliton model (CDM) is Lorentz- and chirally- 
invariant. It has been demonstrated to exhibit dynamical chiral symmetry 
breaking and spatial confinement in the locally uniform approximation. We 
here study the full nonlocal quark self energy in a color-dielectric medium 
modeled by a two parameter Fermi function. Here color confinement is man- 
ifest. The self energy thus obtained is used to calculate quark wave functions 
in the medium which, in turn, are used to calculate the nucleon and pion 
masses in the one gluon exchange approximation. The nucleon mass is fixed 
to its empirical value using scaling arguments; the pion mass (for massless 
current quarks) turns out to be small but non-zero, depending on the model 
parameters. 
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I. INTRODUCTION 



The chromo-dielectric soliton model (CDM) []J is a Lorentz- and chirally invariant low- 
energy effective field theory based on quantum chromodynamics (QCD). In order to simulate 
gluon condensates and other scalar structures (as, e.g., qq pairs) inside hadrons the QCD 
lagrangian density is supplemented by a scalar field a that mediates the gluons through a 
color-dielectric function. 

Following arguments first given by T.D. Lee, a suitably modeled color-dielectric function 
K,(a) guarantees absolute color confinement 0. The assumed potential of the scalar field is 
quartic and has two minima, one at zero and a second, deeper minimum at a finite value 
identified as the vacuum value o~ v . In the absence of quarks, the normal state of the a-field 
is at the vacuum value. In the presence of quarks and gluons, the cr-field finds a minimum 
in the vicinity of zero; the quarks and gluons dig a hole in the vacuum. This is the origin of 
confinement in the model. 

The CDM differs from the original Friedberg-Lee (FL) nontopological soliton model || 
in the essential feature that there is no direct quark-sigma coupling term. Thus the model 
is chirally invariant for massless quarks. Krein et al. || showed that for a locally uniform 
dielectric medium, chiral symmetry is dynamically broken if the strong coupling constant or 
the inverse of the color-dielectric function exceeds a critical value. Consequently, the quarks 
acquire an effective ("constituent") mass. The Nambu-Goldstone boson corresponding to 
this symmetry breaking has been identified with the pion ||. 

While the locally uniform model demonstrated spatial confinement and the emergence 
of the pion, it did not demonstrate color confinement. Furthermore, it was shown that the 
range of non-locality of the quark self-energy was of the order of the typical hadronic length 
scale and hence large compared with the soliton surface. Therefore, it was deemed essential 
to investigate the nonlocal quark self-energy for a realistic and self-consistent soliton. 

This is the problem we address in the present paper. We first obtain the linearized 
(Abelian) gluon propagator in an inhomogeneous color-dielectric medium. Because of the 
Abelian approximation, the calculation is analogous to a problem in electrodynamics. The 
Schwinger-Dyson equation for the quark propagator is solved along the imaginary energy 
axis in order to avoid mass poles on the real energy axis. Quark wave functions are obtained 
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by solving the Dirac equation with the self-energy playing the role of a nonlocal scalar 
potential which is analytically continued to the real energy axis. 

The mutual interaction between quarks and - in the case of mesons - antiquarks in 
hadrons is treated in the one gluon exchange approximation (OGE). Corrections due to 
center of mass motion are taken into account approximately. Using scaling arguments, we 
fix the nucleon mass to its empirical value and calculate the pion mass as a function of 
phenomenological parameters. In the case of massless current quarks, the pion mass turns 
out to be small but non-zero. Since a vanishing pion mass is demanded by Goldstone's 
theorem, the calculated pion mass can be considered a test of the approximation schemes 
applied 0. 

The paper is organized as follows. Sec. [Tl] introduces the basic features of the chromo- 
dielectric soliton model. After deriving the equations for the gluon-propagator in an in- 
homogeneous medium (Sec. Sec. ^ addresses the formulation of the appropriate 
Schwinger-Dyson equation for the quark self-energy. This self-energy is used in Sec. [V] as an 
effective nonlocal quark potential in order to determine quark wave functions in a bag. Sec. 



VI contains details of the numerical solution of the corresponding equations and presents 



results for the self energy. Sec. |V11| describes the calculation of hadronic properties in the 



OGE- approximation and, finally, Sec. |VIII| sums up the main results and discusses future 
prospectives. 

II. THE MODEL 



The CDM lagrangian density is given by [ffl 

Ccdm = q(h"d, + g s 1 -X a A; Y ~ m f )q - k{?)\ F^F^ + \{d,af - U{a) + C , (1) 



F; v = d,Al - dyAl + g s f ahc A\A c v , (2) 



where the color 577(3) structure constants satisfy A a , A = 2if X c , q are the quark fields, 
A^ are the gluon fields, o is the effective scalar field which determines the effective color- 
dielectric function 1 > n(a) > 0, and £' contains any necessary counter terms, gauge 
fixing term, or ghosts. It is evident that the model is gauge invariant, mj is the quark flavor 



(current-) mass matrix. Throughout this paper we will set the current quark masses equal to 
zero, so that the model is also explicitly chirally invariant. The color-dielectric function k(o~) 
mediates the gluon field and is designed to guarantee color confinement. It has been shown 
|| that the following assumptions must be satisfied: k(0) = 1,k((t v ) = K f (a v ) = k'(0) = 0. 
These constraints are satisfied, e.g., by 

k(o-) = l + 6(x)x n {nx- (n + 1)), n > 2, (3) 

with x = <jj<j v . The vacuum- value of the cr-field is denoted by a v . We choose n = 3 for 
simplicity so that n(a) is continuous at x — 0. 

Analogous to the FL-model, the potential of the cr-field is given by the quartic form 

The "bag constant" B is fixed in terms of the other mo del- const ants so that U(cr v ) = 0. 
In the FL-model, U(a) is chosen to be quartic in order to make the model renormalizable. 
Although our model is not renormalizable (due to the presence of ft(cr)) we stick to this form 
in order to minimize the numbers of free parameters in the model. We identify U"(a v ) = vn? GB 
with the lowest ++ glueball mass and require U'(a v ) = 0. 

We now discuss the divergences of the model in more detail. The model exhibits both 
infrared and ultraviolet divergences. The origin of the infrared divergence is the same as 
for the MIT bag model. For a spherical bag, for example, the electric monopole term 
of the quark self-energy diverges as r — > oo. This happens in the CDM, if the color- 
dielectric constant vanishes as r — ► oo. The infrared divergence is thus associated with color 
confinement. However, for a color singlet bag, no infrared divergence occurs since the self 
and mutual interaction terms cancel when ladder diagrams for the mutual interaction are 
properly calculated. The monopole term of the self energy is ignored in the MIT bag model. 
Since this term is the source of color confinement in our model we cannot neglect it. In 
our calculations we choose a Fermi function shaped spherically symmetric color-dielectric 
function as displayed in figure 1. 

The ultraviolet divergence is associated with the point nature of quarks. It is shown in ref. 
Q that the effective quark mass, which is generated because of dynamical chiral symmetry 
breaking, goes to infinity if the color-dielectric function approaches zero. This divergence is 
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thus connected with spatial confinement. We handle this divergence by introducing an energy 
cutoff (asymptotic freedom). For numerical reasons, we regulate the infrared divergence by 
adjusting k v = n{a v ) to a small, non-zero value, and discuss the limit k v — > 0. 



We assume that parts of the non-Abelian effects are effectively included in the a-field. 
This allows us to approximate the gluon field by its Abelian part. Hence the gluon field 
equations are formally identical to Maxwell's equations in an inhomogeneous medium char- 
acterized by a time-independent color-dielectric function k(t). The field equations for the 
vector potential A^(r , t) read (we follow here references and |§): 



Since the Abelian approximation destroys gauge invariance, the choice of gauge is part of 
the approximations. We choose the Coulomb gauge defined by 



III. THE GLUON PROPAGATOR 



&>K(r)[d ll A v -d v A ll ) = J v (r,t). 



(5) 



V- (kA) = 0. 



(6) 



The v = component of Eq. (|5]) yields 



V/e(r) -VAo(r, t) = J {r , t) . 



(7) 



The time-time component of the Green's function, D 00 , defined by 




(8) 



satisfies the equation 



V«(r) ■VD m (r. r') = <f 




(9) 



Note that D 00 (r , r') is instantaneous. 

Now consider the v — i components of Eq. ([5]) 



HiJfA - V 2 (kA ) + V x (kA x VIhk) = J - kV('),A = J 



(10) 



The transverse current defined by Eq. (|T0D can be expressed in terms of J using the time-time 



Green's function: 
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J /r (r./) = J(r./)-K(r)V J d 3 r'D 00 (r, v')d t J (r', t) . (11) 

Using current conservation, d t ,J {] + V- J = 0, and performing a partial integration, we obtain 

J tr (r , t) = J (r , t) - «(r )V J dV (v'D m (r , r ')) • J (r ', t) . (12) 

We now Fourier transform the time dependence of J (r , t) and A (r , t) to J (r , u) and 
A (r , uj). The Green's function corresponding to Eq. ( |T0"D satisfies 

- V 2 + cu 2 +Vx (VhiK) x] K/)(r. r') = ^(r, r') , (13) 

where the components of the transverse delta function are given by 

5%(r, r') = 5 ij 5 3 (r - r') - k(t )d i d' j D 00 (r , r') . (14) 

In this paper we will restrict ourselves to spherical bags. In this case the Green's functions 
can be decomposed in terms of spherical harmonics 0: 

D\, - S(r,r';w) = £ d iB ,(r, r'; ^^(n)?^^) , (15) 
^ 00 (r,r') = Ed?(r,r')^(^;(fl'). (16) 

Imi 

Some of the tensor components are shown in figs. 2 and 3. The n(r) parameters are again 
R = 0.8 fm, A = 0.1 fm, and k v = 0.1. 

It should be noted that the Green's functions do not carry any color indices. This results 
from the fact that the medium is color-neutral so that D^ u has a trivial (diagonal) color 
structure. 

Details of the derivation and solution of equations ([]) and (|IB|) are given in [7J, an 
important correction is reported in ref. ||. 



IV. THE SCHWINGER-DYSON EQUATION IN THE QUARK-GLUON SECTOR 

Being now in the possession of the gluon propagator in the cavity we can study the 
Schwinger-Dyson equation for a single quark in a cavity. In the course of this calculation 
we need not refer to the a-field. 
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The Schwinger-Dyson equation reads (in (u, r )-space): 

E(r , r ; u) = ia' / du'D^r , r ; t/) 7 ^(r , r ; u - u')>f , (17) 

J — oo 

with a' = (4/3)<7g/27r. In Eq. flTTD we have already approximated the one-particle irreducible 
quark-gluon vertex T M by the bare one. 

It is easy to show that both the gluon propagator D^ v and the quark propagators G do 
not have poles off the real cu-axis. So according to the Schwinger-Dyson equation, the self 
energy E should have no pole of the real a;-axis as well. Thus we perform a Wick rotation, 
oj — > iy and study the self energy first for imaginary u. The "rotated" Schwinger-Dyson 
equation now reads: 

E(r , r ; y) = -a' / dy'D^r , r ; tf)<fS(v , r ; y - y')Y ■ (18) 

Simultaneously, the Dirac equation for the quark propagator has to be satisfied: 

(^7°-7-p-E)S = 5 3 (r -r'). (19) 

To simplify the notation we have used the shorthand US for / (i 3 r2E(r , r 2 ; uj)S{t 2, r '; cu). 
Throughout this paper repeated spatial coordinates are integrated over. 
We now define the following hermitian functions, 

G = -S(3, V = PE . (20) 

Eqs. ( |T9"D and ( TTT] ) then become 

(-uj +a-p + V)G = 5 3 (r -r'), (21) 

V(r , r ; y) = a' / dy'D^v , r ; y')a»G{v , r ; y - y')a v . (22) 

J — 00 

Here a M = (l,a) is used for notational convenience only. It is obviously not a Lorentz vector. 
From the coupled Eqs. (21, ^2|), the hermiticity of G and V can be verified: 

Gt(r,r» = G(r,r';^), (23) 

V i (r,T';u) = V(T,T , ;u*). (24) 
The hermitian conjugation includes the interchange of the arguments r and r': 

Vj.(r,r';a;)=^(r',r; W r. (25) 
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A. Angular decomposition of the quark propagator 



In order to solve the coupled Eqs. ([H], |22|) numerically, we make an angular decomposi- 
tion of the appropriate quantities. For spherically symmetric color-dielectric functions k(t), 
the hermitian properties G and V can be decomposed P,|10| 



G(r,r'; W ) = £ 



g^{r,r , \u))'K K g 12 (r,r';uj)ia r iT_ K 
-gf(r,r';u)ia r 7r K # 22 (r, r'; w)tt_ k 



(26) 



V{v,v'-u) = Y. 

K 

where the respective angular part is given by the 2x2 matrices 



k \ — f 21 (r, r'\ u)ia r 7r K t> 22 (r, r'; c<j)7r_ 



(27) 



The following reduction relationship holds: 



(28) 



(29) 



3^ KAt (f2) are the usual two component spinor spherical harmonics. They are eigenstates of 
the operators J 2 , L 2 , J 2 , and # = ( J + l/2)(-l)( J - L+1 / 2 ) [| 



m;,m s ^ 



(30) 



and obey the orthonormality relation 



Jdnyl(n)y K ,^(n) = 5 KK ,5^. 



(31) 



The radial functions g and v have the following symmetry properties 



g^r, r'; u) = g^r', r; d) = g»(r, r'; u;*)* . 



(32) 



Inserting Eqs. (|26|) and (p7|) in the Dirac Eq. for the quark propagator and using 



(0) and (Ea) yields 



— 



— 1 / r — d/dr + n/r 



1/r + d/dr + n/r 



+ 



^ n ^ 12 
^ 21 ^ 22 
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5(r — r' ) 



rr 



(33) 



where vg denotes / r 2 dr2v(r,r2]u>)g(r2,r';uj) for notational convenience. 

Defining g(r, r'; us) = rr'g(r, r'\ u) and v(r, r'\ lS) = rr'v(r, r'; u) Eq. ( |33"D simplifies finally 

to 

— —d/dr + k/t 



vi 1 e 



d/dr + n/r —uj ) \v 2 } v 22 



K K 



gl 1 gf 
e gl 2 



5(r - r' ) . (34) 



Details of the non-trivial numerical solutions of this equation are discussed in Sec. VI 



B. Radial part of the Schwinger-Dyson equation 



Inserting (|26| ) and ( p7|) into Eq. (0), we can write 



V(r,r ; y) 



v}. 2 (r, r'\ uS)i(j r ix_ 



-vl 1 (r,r';u)ia r TT K v 22 (r,r';uj)n_ K 



a' J dy'D^r , r '; y')a"G{r ,r';y- y')a l 



(35) 



,. ^2(r,r»7r_ K -^i(r,r';^a r 7r 

-a \ dyd jlv {r,r ■y)a-y jlm \ „, , ^ 3^™'°" 

, £ K (r, r ; u;)w7 r 7r_ K # K (r, r ; ^)vr K 



+ a' j dy'd^ry^y'^ 
From the Appendix we find 



^(ry^vr,, g 12 (r,r';u)ia r n_ K 



v 
1 1, 



g K {r,r ;uj)ia r 7T K g z K z {r,r ;u))it- K 



(36) 



(37) 



(38) 
(39) 



The self-energy coefficients A, B and C are explicitly defined in the Appendix. 
With these formulae the radial Schwinger-Dyson Eq. (|36|) reads 

v l R \ry-y) = a' [ dy'\d°(r, r')^!(r, r'; y')Cf 
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- d jW (r, r'; y')g%{r, r'; y - y')A%, K ' 



(40) 



v K (r, r'; y) = a' dy' 



d°(r,r'W(r,r r ;y')Cr K - K ' 



+ dMr,r';y')g 2 J(r,r';y-y')B- i r 



(41) 



v K (r, r'; y) = a' dy' 



d°(r,r')gl}(r,rW)Cr' 



+ d w (r,r / ] y')g 1 J(r,r' ] y-y')B^' 



(42) 



vl 2 (r,r';y) = a' dy' 



d°(r } r')fJ(r } r';y')Cr- K ' 



d jlv (r,r';y')gl}(r,r';y - y')A^ 



(43) 



Here we have also used Eqs. flTo] ) and flToT) as well as the symmetry properties d^{r, r') = 
d®(r',r) and djw(r,r';uj) = dji>i(r',r;uj) which hold for both pure real and imaginary uo 
because in Eq. (|l"3"D for the gluon propagator only u 2 (and not a;) appears. 

V. THE QUARK WAVE FUNCTION 



By interpreting the nonlocal self-energy as an effective potential we can now determine 
the wave function q(r) and energy eigenvalue e of a single quark in the cavity. The corre- 
sponding Dirac equation reads 



a p q(r) + / d 3 r 2 V(r , r 2 ; e)q(r 2 ) = e q(r) . 



In spherical coordinates, q(r) can be written in the form 

u K (r)/r 



q(r) = ^2 

w \-ia r v K (r)/r 



After angular decomposition, the radial part of Eq 

—d/dr + K,/r\ ( u K (r) 

d/dr + n/r / \v K (r) 



obeys 



dr 2 V K (r, r 2 ;e) 



(44) 



(45) 



u K {r 2 ) \ ( u K (r) 
v K (r 2 ) J \v K (r) 



• (46) 
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VI. NUMERICAL CALCULATION 



In our calculations we use a (modified) Fermi function shaped spherically symmetric 
color-dielectric function 

1 — K 

K ^ = l + e (r-R)/A + > ( 4? ) 

where R and A are the radius and the surface thickness respectively of the profile (see fig. 
1). The small but non-zero vacuum value k v guarantees that, e.g., the energy of a single 
quark in the cavity remains finite. For color-singlet multi-quark systems the limit k v — > 



can be performed as will be shown in Sec. VTJ 



Since our model is not renormalizable, an ultraviolet momentum cutoff is needed; this 
is consistent with asymptotic freedom. This cutoff should reflect the energy scale of the 
described physics; we choose Aqdm = 5.0 fm _1 . In terms of the variables used in this paper 
the u/-integration in Eq. ( |P7|) is cutoff at |o; maa; | = Acdm an d the necessary summations over 



angular momenta are limited by l max = R ^max with R from Eq. (|47|). A careful analysis 
of the renormalization problem for a nonlocal, spatially varying dielectric medium can be 



found in ref. 11 



With the gluon propagator derived in Sec. [TTT], we solve the coupled equations ( |3"4j ) and 



( flOl - ^) to obtain the full quark propagator and the quark self-energy on the imaginary 
cu-axis. Because of the absence of poles in this region the self energy is numerically stable 
and no oscillations occur. A Taylor-expansion method is subsequently applied to construct 
the quark self-energy on the real u>-axis: 

^2 Z^ 

v K {r, r'; z) = v K (r, r'; 0) + zv' K (r, r'; 0) + — <(r, r'; 0) + — (r, r'; 0) + — (r, r'; 0) + ■ • • , 

I o 24 

(48) 

where the derivatives are evaluated in terms of the discrete values of the functions along the 
imaginary u;-axis. 

Eq. ( 151 ) is a coupled system of integro-differential equations. For its numerical solution 
we use matrix inversion. It is well known that the leap frog instability [|l^] (p. 342) appears in 
an equation like Eq. (|34]) when the first order derivative is replaced by a centered difference. 
Therefore we introduce a small second order derivative term to suppress the instability: 
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-OJ 



-d/dr + Bd 2 /dr 2 + n/r 



d/dr + Bd 2 /dr 2 + n/r 



-co 









)* 







a 11 a 12 

n 21 n 22 
, UBk iJBK , 



5(r — r' 



(49) 



where B = b A is a small number. A is the grid interval and b ~ ±1 for k — ^1 (the sign 
is important to suppress the leap frog effect!). This additional regularizing term does not 
spoil the accuracy of the solution. 

Numerically we find that g K = g% K satisfies Eq. (|33|) very well if B ~ ±A. Therefore 
qbk, can be considered as a first approximation to g K . There must be a discontinuity in the 
Green's function for first order differential equations. In Eq. ([34]), this discontinuity occurs 
in its off diagonal elements [|nj. We find numerically that the additional second derivative 
term smooths ot the discontinuity somewhat. 

This approximation can be improved. This will be demonstrated first in general terms. 
Consider the following two Green's equations: 



(L + L B )G = 6. 

After operating with Go on both sides of Eq. (^l|) and integrating, we have 

Gq = Gq(L + Lb)G 
= G + GqLbG , 

or 



(50) 
(51) 



(52) 



G — Go — GqLbG . (53) 
Similarly, by integrating both sides of Eq. ([£|) with we have the exact relation: 



g K (r,r';uj) = g BK (r,r';uj) 



+ J dr 2 g 0K {r, r 2 ;uj) 



Bd 2 /dr 2 



Bd 2 /dr 2 
12 



v K (r 3 ,r 2 ;u})g K (r 2 ,r';u) . 



(54) 



This equation can be solved by iteration. However, in this case the leap frog instability 
eventually creeps in again. We have thus carried out only one iteration. 

Some of the results of the self energy calculation is shown in figs. 4, 5 and 6. The figures 
display V™ n (u;r,r') for k = — 1, u = lfm" 1 and (run) = (11), (22) and (12) respectively. 

We note that the self energy is non-zero. This implies a dynamical breaking of chiral 
symmetry. The structure of the self energy reflects the nonlocal character of the interaction. 
However, the self energy is sharply peaked around r — r' reflecting the dominance of the 
local contribution. 

The self energy is inserted in the Dirac equation which is solved self-consistently for the 
ground state (k = —1). The result is shown in fig. 7 for the n(r) profile of fig. 1. 

The single-quark energy e is shown in fig. 8 as a function of k v . It does not exhibit a 
sign of divergence as far as the calculation could be carried out (down to k v = 0.05). In 
fact, e turns out to be quite insensitive to k v if k v is small. The presented results are thus 
gratifying. 

We have tested our numerical calculation by varying the following numerical parameters: 

(a) the number N max of r grid points, 

(b) the integral limit r max of r and 

(c) the number N u of uo grid points. 

For the actual parameters we have chosen the physical observables are all insensitive to 
them. 

VII. HADRONIC PROPERTIES 

Having calculated the wave function and energy eigenvalue of a single quark in a cavity, 
we now investigate color-neutral composite systems of N q valence quarks. Evidently, N q = 2 
for mesons and N q — 3 for baryons. The energy of these systems is calculated in the one 
gluon exchange approximation. Finally, corrections due to the center-of-mass motion and 
to the a-field are taken into account approximately. Using scaling relations we fix the mass 
of the nucleon and study the pion mass m n . Its deviation from zero is a measure of how 
good our approximations are since the pion should be massless according to Goldstone's 
theorem. 
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A. One gluon exchange approximation 



The one gluon exchange interaction energy between quarks (of equal eigen-energy) is 
given by @ 

E ex = a' J d 3 ri d 3 r 2 [f(r 1 )D 00 (r h r 2 )j°(r 2 ) - j (r x ) • £(r i, r 2 ; 0) • j (r 2 )] , (55) 

with a' = \gl J2i<j <\ ' \ >■ The color matrix element <Ai • A 2 > has the value —16/3 for 
the pion and —8/3 for the nucleon [0. Taking into account that that for both the pion and 
the nucleon, the quarks are in the ground state with n = = ±1/2 the corresponding 

currents can be evaluated and the exchange energy is ready calculated. 

The total energy of quarks and gluons in a hadron with N q valence quarks is then given 

by 

E q>g = N q e + E ex . (56) 



B. Corrections and sigma contributions 

Up to now the a-field has been neglected. However, it contributes to the total energy of 
the bag. The cx-field can be reconstructed from n(r) and n(a) given in Eqs. (f47D and (^|) 
respectively. Then the cr-field energy is given by 



E n 



U(a) 



(57) 



with the potential U(a) given in Eq. (^). The total energy of the bag is then Eb ag = E q ^+E a . 

We now address the hadronic center-of-mass energy. Since localization of the bag breaks 
Lorentz invariance, the bag acquires a non-zero total momentum that contributes to the 
total energy of the system. The easiest way to correct this effect is to use the following 



approximate formula (projection |T3[ would be better but more cumbersome): 



m. 



E 



bag' 



<P 2 > 



bag j 



(58) 



with 



< P 2 > 



bag z 



N q < P 2 > a + < P 2 > CT . 



(59) 
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The momentum squared of one quark is given by 

< p 2 >q = J d 3 r\Vq\ 2 . (60) 

In order to calculate < P 2 > a , the coherent state approximation || is used: 

< P 2 > a = J d?kk 2 u k f 2 . (61) 

fk are the Fourier transforms of er(r), u>k is the cx-field energy in the mode k. 
For slowly varying Uk we finally get: 

< P 2 > a = (m 2 GB + < k 2 >)5 J d 3 r(Va) 2 , (62) 

with the glueball mass uigb- 

C. Scaling and the nucleon mass 



Scaling can be used to generate new solutions [[L4J] from those presented so far. The 



equations are invariant under scale transformations where all lengths r are replaced by 

r — > r' — Xr , (63) 
all energies and frequencies (including the cutoff Acdm) are replaced by 

E — > E' = E/X (64) 

and 

a a' = a/A 2 , 6' ->6/A. (65) 
c and a s are invariant. The a-field and the gluon field potentials scale as length -1 . 

D. Numerical results 

Throughout our calculations, we use a cutoff Ac dm = = 5 fm -1 . With l m — R u m , 
the quark wave functions and energies depend on the two parameters R and A from the 
k(t) profile. The hadron masses additionally depend on the parameters a, b, c of the a-field 
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potential U{a). In order to minimize the number of free parameters, we assume that U(cr) 
is universal in all hadrons. However, each hadron has a different /t(r) profile reflecting the 
fact that the hadronic size is not universal. 

The numerical procedure is as follows: We choose a potential U(a) and calculate the 
corrected nucleon mass according to Eq. (|5Bp. Using scaling relations we renormalize all 
dimensional properties by fixing the nucleon mass to its empirical value = 938 MeV. 
With these renormalized properties we now calculate the pion mass as a function of R and 
.4. 

To this point, the cr-field is not self-consistent. We now vary the parameters of the n(r) 
profile in order to find an extremum in the energy. This is a first approximation to a fully 
self-consistent treatment. However, we expect the results to be reasonable since the proper 
shape of the cr-field is similar to a Fermi function. 

We find that there is not always a minimum in A for a given m n (R). This may be due 
to the crude method used to correct the effects of the-center-of mass motion, to the form of 
ft(r), or the point is an extremum, not a minimum. The resulting pion mass is small but 
non-zero. 

VIII. SUMMARY AND PROSPECTIVES 

Within the framework of the chirally-invariant chromo-dielectric soliton model, the 
Abelian gluon propagator is solved in configuration space for a color-dielectric function with 
two parameters. The quark self energy was obtained by solving the (nonlocal) Schwinger- 
Dyson equation in configuration space as a function of imaginary energy. Quark wave 
functions and real eigenvalues were obtained. Bag states were constructed for the pion and 
the nucleon including one gluon exchange mutual interaction between quark pairs. The 
parameters of the parameterized cr-field (or equivalently, the dielectric function n(r)) were 
varied to extremize the bag energy. Approximate center-of-mass corrections are calculated. 
Employing scaling relations, the nucleon mass was set to its empirical value. The resulting 
pion mass was determined to be small (the actual value depending on the model parameters) 
but not zero, as demanded by Goldstone's theorem. 

Extensions of the present work include the following: 
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a) Center-of-mass corrections based on variation after projection. This technique has 
been studied extensively by Liibeck et al. []I3| for the Friedberg-Lee soliton and was found 



to give significant corrections. It is certainly more reliable than the prescription m? k-< 
H > 2 — < p 2 > used in the present paper. 

b) A "more" self-consistent treatment of the soliton field by either solving the differential 
equation for the a-field or by including more parameters in the functional form of n(a) and 
k(t). 

c) Calculation of the mutual gluon exchange between quark pairs by full summation of 
ladder diagrams. 

d) A systematic adjustment of model-parameters to fit the properties of all low-lying 
hadrons. This is not as tedious a task as it might first appear. The parameters of the model 
are a, b, c and a s . The functional form of n(a) also introduces a model dependence, but the 
results appear to be quite insensitive to that. Of the four, one is set by the nucleon mass 
using scaling from any given set. Results appear to be relatively insensitive to the "family" 
characterized by b 2 /ac @ but this is related to the glueball mass which is assumed to lie in 
the range of 1 — 2 GeV. The key parameters include the nucleon size, magnetic moments, 
9a/ 9v an d the N-A mass splitting. Other hadronic spectra properties are then regarded as 
predictions of the model. 
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APPENDIX A: 



In this appendix the self-energy coefficients A, B and C from Sec. [IV B| (Eqs. fl37r[39|) ) 
are explicitly evaluated. 

We start with formula (5.9.15) of Edmonds' jnj 



o q Xv = VS< 1/2, u, 1, q\l/2, q + v> X q -\ 



(Al) 
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Note that we use throughout our calculations the phase convention of Edmonds. With the 
definitions 



(A2) 



'm'£q , 



<jrm' 



(A3) 



x ± iy 

e ±i = =F ^ ' c o = z > 

we get 

= E <l '' m i> !> 171 > Y i'mi<r q < m 2 , 1/2, v\j K , fj, > XvYl, 



(A4) 



m,2 



E < I', mi, 1, m > ^lm < "^2, 1/2, Z 1 >< I', 0, Z K , 0|L, > 

miqm2vLMv' 



x < mi, l K , m 2 \L, M > 



/ 2 ™i| 1) ^ <lAM ' 9|1/2 '^ > ^ 



= E V( 2/ + + 1)(2/ K + 1)(2Z' + 1)3/2tt 

' L 1/2 y ^ 

x <l,m, j K , jj! >< I', 0, l K ,0\L,0 > < 



> y , 



(A5) 



V 1 Z 
U 1/2 j K ) 

Here we have used the contraction formula for spherical harmonics (Edmonds (5.16)) and 
the definition of the 9j-symbols (Edmonds (6.4.3)). 
Similarly, 



/47T, 



E 2 < I', mi, l,q\l,m > Y Vmi \ —Y lq < l K , m 2 , 1/2, u\j K , fi > XvYi K m,_ 



= E 2 <r, 0, 1,0|Z,0> Y lmt 



m\qm2V 



4tt (2 + l)(2Z'+l) 
\ 3 4tt(2Z + 1) 



< l K , m 2 , 1/2, Z/|j«, /i > 



K,m 2 



E( _ 1)1/2+ , + , +j /(2i. + l)(2^ + l)( 2y + l) < r> _ t _ 0|; _ 



> 



x < Z,0,Z K ,0|L,0 > 



j« 1/2 l K 
L I j 



<i,m, j K , fi\j, fi' > y^, 



(A6) 
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and 



E Y lm < L, m 2 , 1/2, u\j K , \L > XuYl K m 2 



'(2L + l)(2j K + l)(2/ + l) 
Arc 



x < Z,0,Z K ,0|L,0 > 



< 3k, /J, i, m\j, n' > y^> ■ 



(A7) 



j K 1/2 Z« 
L I j 

According to Eq. ([AJ), the expression E^DW^) • a3^(ft)3^(ft> ■ 3?,„ m (fi') is 
proportional to 5jj> and Now L and L' have to be equal or differ by 1. However, 

< I', 0, l K , 0\L, >< Z", 0, Z K ,0|L', > vanishes if |L — L'\ is odd, since I' — I" is even, so 
only terms with L = L' (or k' = k") contribute. Thus 



E%'m( n )-^(n)3^(n'>r.3>£ Wm (n / ) 



(A8) 



The following symmetry relation Afi/f/i = Af^, holds. 

Similarly, according to Eq. (gg), the expression £ mM 2f ■ 3^/ m (f2)34 M (fi)3^(fi')<x • 
3^// m (fi') is proportional to <5.,y and Now L and L' have again to be equal or dif- 

fer by 1. However, < Z, 0, l K , 0\L, >< V , 0, 1, 0|Z, >< I", 0, l K , 0\L', > vanishes if \L-L'\ 
is even, since I' — I" is even, so only terms with L = L' ± 1 (or k' = —k") contribute. Thus 



E4n"^ v (^)^v(^)- 



(A9) 



K/t' 



Similarly, we have 

E^m^(^)^L(^)^(^) = E^^V^^v^') 



(A10) 



Finally, the quark-gluon coupling coefficients ^4, £> and C are given by 
3 



•A-ll'l" 



2vr 



(2/' + 1)(2Z" + l)(2j« + 1)(2Z„ + 1)(2Z + 1) < 0, l K , 0\l K ,, > 





'In' 1/2 




1/2 




x < l",Q, Z K ,0|Z K ,,0 > < 


Z' 1 / 


II 


Z" 1 Z 


(AH) 




Jk 1/2 j«. 




A 1/2 J K , 
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Y3K K 

°ll'l" 



■1) 



-l/2+l+l K +j K , 



3/2(2/' + 1) (21" + 1) (21 + 1) 



(2/ K + l)(2j K + l) 



7T 



x < Z", 0, l K , 0\l K ,, >< I', 0, 1, 0|Z, >< 1, 0, Z re , 0|Z K -,, > 
f Z re / 1/2 j K > ) 

3k 1/2 Z re ) 

U 1/2 J J 



(A12) 



cr ^(2^ + l)(2j. + l)(2Z + l) < 2 fi. 1/2 l H 

47r / j K > 



Furthermore, 



E[2r .^(n) - a^ ra (n) -o]^(n)yt (n>. ^(jy) 



my, 



In that very last step we have used Eqs. 
and the identification 



and 



(A13) 



(A14) 

as well as the identity a T y Kfl = —y R y 



Duty i = Duly i 



ll'l" ~ ^ll'l" 



(A15) 



Working out the hermitian conjugate of Eqs. ( |A"8| ) and ( |A^| ) we get the following symmetry 



relations: Afyf,, = A%% and B%% = B„,?„ K 



Hl"l 



wi'' 



J ll"l' 
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Figures 



Figure 1: The color-dielectric function k(t) for R = .8 fm, A = .15 fm and k v = .15 (r 
in fm). 

Figure 2: The tensor part of the gluon propagator in the transverse magnetic mode 
dio2(r,r'). 

Figure 3: The tensor part of the gluon propagator in the transverse magnetic mode 
di22(r,r'). 

Figure 4: The quark self energy on the real a>axis: v}\ (r, r') . 
Figure 5: The quark self energy on the real tu-axis: v 2 \ (r, r') . 
Figure 6: The quark self energy on the real u;-axis: v}\ (r, r') . 

Figure 7: The quark wave function. ru(r) is the darker line, rv(r) is the lighter one. 
Figure 8: The single-quark energy e (in fm -1 ) as a function of k v 
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TABLES 

TABLE I. Table of pion masses for a = 39.9, b = -746.2, c = 4569.6, B = 0.03892, 
m GB = 1310.8 

















R 


A 


E q 


E q,9 


V< P 2 >Q 


V< P 2 > CT 


m w 


(fm) 


(fm) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


0.6 


0.150 


404.92 


610.61 


477.20 


189.6 


171.31 


0.6 


0.175 


405.81 


612.17 


476.50 


178.3 


235.94 


0.6 


0.200 


405.97 


613.19 


476.10 


170.1 


290.23 


0.6 


0.225 


405.71 


615.40 


475.90 


164.0 


344.06 


0.6 


0.250 


405.22 


614.81 


475.80 


159.4 


390.66 


0.8 


0.150 


295.37 


452.89 


489.10 


249.4 


344.64 


0.8 


0.175 


296.77 


454.19 


486.50 


232.1 


284.70 


0.8 


0.200 


297.20 


454.17 


484.80 


219.0 


208.37 


0.8 


0.225 


297.01 


453.40 


483.70 


208.8 


61.61 


0.8 


0.250 


296.39 


452.22 


482.90 


200.7 


198.74 


1.0 


0.150 


230.69 


358.91 


485.20 


311.7 


338.67 


1.0 


0.175 


232.25 


360.33 


480.00 


288.5 


238.99 


1.0 


0.200 


232.76 


360.09 


476.20 


270.6 


27.71 


1.0 


0.225 


232.60 


358.83 


473.80 


256.3 


237.19 


1.0 


0.250 


231.99 


356.95 


471.70 


244.8 


342.74 


1.2 


0.150 


194.97 


308.29 


400.90 


375.5 


371.10 


1.2 


0.175 


196.45 


309.63 


398.10 


346.6 


462.86 


1.2 


0.200 


197.08 


309.56 


396.10 


323.9 


540.69 


1.2 


0.225 


197.18 


308.61 


394.60 


305.7 


612.00 



23 



